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We combine model mapping, exact spectral bounds, and a quantum Monte Carlo method to 
study the ground state phases of a mixture of ultracold bosons and spin-polarized fermions in a 
one-dimensional optical lattice. The exact boundary of the boson-demixing transition is obtained 
from the Bethe Ansatz solution of the standard Hubbard model. We prove that along a symmetry 
plane in the parameter space, the boson- fermion mixed phase is stable at all densities. This phase is 
a two-component Luttinger liquid for weak couplings or for incommensurate total density, otherwise 
it has a charge gap but retains a gapless mode of mixture composition fluctuations. The static 
density correlations are studied in these two limits and shown to have markedly different features. 



Recent advances in experiments with cold atoms in 
magneto-optical traps have led to a series of exciting 
results. One area of emerging interest is the study of 
degenerate mixtures of bosons and fermions in various 
trap and lattice geometries. While the original motiva- 
tion for such studies was the sympathetic cooling of the 
fermions via interactions with bosons, such systems pro- 
vide a unique opportunity to study experimentally quan- 
tum phases and associated transitions in a system of 
mixed quantum statistics, which are rare in condensed 
matter systems. The absence of defects and impuri- 
ties as well as the ability to tune the interactions be- 
tween the particles will enable the exploration of various 
quantum many-body phenomena. Stable boson-fcrmion 
mixtures have been realized with 6 Li- 7 Li|l|], 23 Na- 6 Li[^), 
87 Rb- 40 K@ and 6 Li- 87 Rbf§. On the theoretical front, 
the ground state properties of boson-fermion mixtures 

in the continuum and on a lattice - have been in- 
vestigated using mean field theories 0, ||, |7|, |[ |9[ [l0| , 
bosonization|llJ, Bethe ansatz (BA)pH and numerical 
simulation of small systems (lj, H . 

In this work we study the nature of the ground state 
of a boson-fcrmion mixture on a ID optical lattice, con- 
centrating specifically on the special case of equal hop- 
ping coefficients for bosons and fermions, ts = tp in 
Eq. (0). We construct mappings to several known mod- 
els in appropriate limits and use known analytic re- 
sults from these to supplement our numerical results. 
In addition to exploring the unique quantum demix- 
ing transitions, we have investigated the nature of the 
ground state in the homogeneous (mixed) phase at dif- 
ferent interaction strengths and densities. Our results 
are summarized in Fig. |l|. At large fj,p (^b), the ground 
state is purely fermionic (bosonic) with nsinp) = 0, 
whereas a thcrmodynamically stable mixed phase (B-F) 
with ub, tif > exists over a finite range of parameters 
around fj.s ~ The exact analytic boundary between 
the pure fermion F-F and B-F phases obtained from the 
solution of Eqs. (|J) and (|J) is also shown. The agreement 
between the two approaches is excellent in the weak cou- 
pling limit, whereas the B-F phase is overestimated in 
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FIG. f: Ground state phase diagram of the boson-fermion 
model ([j]). Phases are purely bosonic (B-B), purely fermionic 
(F-F), and a uniform mixture of the two (B-F). The filled 
squares joined by dotted lines to guide the eye represent the 
QMC results, while the solid lines denote the analytic phase 
boundaries. The dashed lines separate the commensurate 
(charge-gapped) and non-commensurate (gapless) sectors of 
the corresponding phases. In the weak coupling limit, the 
B-F and B-B phases do not have a charge gap. 

the numerical data at strong coupling. The extent of the 
mixed phase gets smaller at low densities, especially at 
large U and V. For the case U = V, we show that even 
at very strong coupling, there is no direct transition be- 
tween pure bosonic and fermionic phases. A mixed phase 
is always stable, contrary to the mean field results 0. At 
weak couplings, the mixed phase is dominant and one 
needs large differences in the chemical potentials for the 
two components to observe demixing, consistent with the 
bosonization results |ll|. 

A mixture of interacting bosons and spinless fermions 
in an optical lattice is well described by the one-band 
Bose-Fermi Hubbard model (BFHM), 
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Here, U and V parametrize the on-site boson-boson and 
boson-fermion interactions; we assume U, V > 0. The 
chemical potentials fig and \ip couple to the densities 
of the individual components, ng and np. In the experi- 
ments, there is also a trapping potential that confines the 
particles to a certain region of the lattice. As a simplifi- 
cation, we have ignored such a confining potential. The 
present description is applicable in the central part of a 
system with a shallow trapping potential. In the discus- 
sion of the numerical results, we assume ts ~ tp = t = 1 
which sets the energy scale. All other parameters are 
(implicitly) expressed in terms of t. 

An insight into the model (Q) can be gained by study- 
ing the limiting cases. For V = 0, the boson and fermion 
sectors are decoupled. The fermions are non-interacting; 
they fill a band with momenta |fc| < kp = imp. The 
fermion sector has a commensurate density and an associ- 
ated band gap for > 2tp; otherwise it is gapless. The 
boson sector is described by the Bose-Hubbard model |l5|. 
At incommensurate densities ns, the ground state is a 
gapless supcrfluid (SF) for any U > 0; at commensurate 
densities, there is an SF to Mott insulator transition at 
U rs 3.5tB [fl6|| . Finally, the boson sector is trivially an 
insulator at zero density, rig = 0, at [Ib < —2ts- 

In the hard-core limit U — > oo, the bosons can be 
mapped to a second "flavor" of fermions. Then, for 
ts = tp, the resulting model can be rewritten as the 
usual Hubbard model, where the fermion flavor serves as 
the (pseudo)spin index. The effective external magnetic 
field is h = fip — fig, the Hubbard coupling constant is 
U = V, the hopping t = ig, and the chemical potential 
ft = (Ms + Mf)/2 couples to the density n = ns + np. 
This model is exactly solvable ]l8||. Most importantly, 
there is no charge or spin gap away from the commensu- 
rate density, < n < 1. Given the density n, the Hub- 
bard fermions reach full spin polarization at the field JTs|] 

7 , U f™ „ , COS k — COS TTTl 

h c (n,U) = — dkcosk—, — z — . (2) 

y ' it J {Jj/Atf + sin 2 k y ' 

At this point one of the Luttinger components disap- 
pears, we are left with only one kind of the fermions. 
The corresponding value of the chemical potentials 

fiF — —2t cos irn, ft = fip + h c (n,U)/2. (3) 

The system acquires a gap at the commensurate density, 
n = 0,1. The corresponding values of the critical mag- 
netic field are h(0, U) = 0, htt, U) = (U 2 + 16?) 1/2 - U. 

For the original BFHM (@), the field h = Th c (n,V) 
would correspond to the demixing transition, where the 
density of either fermions or bosons turns to zero. We 
also note that in the latter case, ns = 0, the boson re- 
pulsion constant U is irrelevant. Therefore, as long as 
the transition is continuous (as expected for V < 2[7), 
the boson demixing transition, ub = 0, is expected at 
jjLB = (J*f ~ h c (n,V). The corresponding curve on the 



^f~^b diagram [Fig. [l| is given parametrically in terms 
of n = n F by Eqs. (|), (|). 

To understand the fermion demixing transition at 
np = 0, notice that in the "symmetric" case, ts = tp, 
U = V, and /is = the Hamiltonian (Q) commutes 
with the superoperator, Q = YlfJbi, an d its conjugate, 

. A ground state wavefunction (WF) with N bosons 
and no fermions, can be mapped to an eigenstate 

with N - 1 bosons and Np = 1, = N^ 2 
The normalization was found using the anticommutator, 
{Q, Q t } = N B + N F . While may not be the ground 
state, the degeneracy implies that adding fermions is def- 
initely favorable for (ip > fiB- Combined with the local 
instability of the fermion phase at the line Eqs. (|2[), (^), 
this gives a proof of thermodynamical stability of the 
mixed phase just above the symmetry line. 

The same argument could also work the other way, as 
long as the ground state with Np = 1 fermion, &bi, 
contains a fully symmetrical component, so that the cor- 
responding bosonic WF Q^bi ^ ®- This happens at 
strong coupling (certainly at n = 1 discussed below), as 
well as at small densities, as can be seen by analyzing the 
energy of the continuum BA solution[|l2| (a). In these two 
cases, the fermion demixing transition is indeed located 
precisely at the line /ip = \xb (solid diagonal in Fig. 1). 

Finally, for strong but non-infinite couplings, U, V 3> 
tB,tF, the low-energy physics of the system is described 
by the effective Hamiltonian which at ts = tp is analo- 
gous to the anisotropic t-J-V or extended Hubbard model 
in external magnetic field [^0[ . The model has an easy- 
axis anisotropy for V > 2U , which implies a spin gap dis- 
appears via a first-order "spin-flop" transition (the dis- 
continuous demixing transition is accompanied by phase 
separation). At the commensurate total filling, n = 1, 
the system acquires a charge gap and is further reduced 
to the exactly-solvable XXZ spin-chain with the param- 
eters J z = 2(t% + t F )/V - t 2 B /U, J x = -2t B t F /V[^ 
in an additional magnetic field so that ns = 1 at the 
symmetry line discussed above. The operator Q under 
this mapping becomes a zero-fc spin-wave creation oper- 
ator. Again, there is no spin gap for |J | > J z (the BF 
mixture is a single-component Luttinger liquid), while 
for | J- 1 1 < J z the system is fully gapped and the corre- 
sponding fermion demixing transition is first-order. 

We have used the Stochastic Series expansion (SSE) 
quantum Monte Carlo (QMC) method|2^] to simulate 
the BFHM in the grand-canonical ensemble on chains of 
length L = 8 — 128. Ground state results are obtained 
by taking large enough values of the inverse temperature, 
(3, where j3 = 2L was sufficient. To characterize different 
phases, we have studied the charge (supcrfluid) stiffness 
of the fermions (bosons), p« = (VF 2 )/2/3L[||, a = F,B, 
and the static susceptibility for the density correlations, 
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FIG. 2: The density of bosons and fermions as a function of 
the bare fermionic chemical potential, /xp. 

We start the discussion of the numerical results by 
identifying the range of parameters which result in a ther- 
modynamically stable homogeneous mixed phase. We 
consider fixed (and equal) values of U and V and vary 
\if and hb- We have analyzed two cases in detail - one 
in the weak coupling regime (U = V = 2), and the other 
in the strong coupling regime (U = V = 10). Since 
SSE is formulated in the grand canonical ensemble, one 
never observes a ground state with spatially separated 
domains. Instead, depending on the parameters of the 
Hamiltonian, the ground state will be in one of three 
phases discussed above - B-B, F-F, or B-F. Results for 
L = 32 — 128 are shown in Fig. ||, where the bosonic and 
fermionic densities are plotted as a function of varying 
fip at a fixed value of (is- The upper (lower) panel shows 
the results for the weak(strong) coupling regime. With 
increasing \ip (at fixed pb), the ground state evolves con- 
tinuously from B-B to B-F and finally to F-F phase. The 
data (along with those for other values of /is) are com- 
bined to obtain the phase diagram Fig. 1. 

Experimentally, the transitions can be observed in a 
large system with fixed particle numbers if the trapping 
potentias seen by the two species of particles vary dif- 
ferently. For example, one can have a shallow trapping 
potential for the bosons and a faster varying one for the 
fermions. In that case, there will be purely fermionic do- 
mains at the edges and a uniform mixture at the center 
of the trap. Our data show the transitions to be continu- 
ous, in contrast to the mean field results, where the B-F 
to F-F transition is found to be discontinuous [0 . 

Next we investigate the properties of the mixed phase 
in detail across the parameter range explored in Fig. [|. 
We start with the charge (phase) stiffness of the fermions 
(bosons) - a finite value implies a gapless spectrum for, 
and consequently, a Luttinger liquid character of, the 
fermionic (bosonic) sector. The results arc shown in 
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FIG. 3: Charge stiffness of the fermions (pf ) and thephase 
stiffness of the bosons {pf ) for the parameters in Fig. 0. 

Fig. |^. At weak couplings, pg {p$) is finite for any non- 
zero us (it-f) (except pf = for tip = 1 as the fermions 
form a filled band). In particular both pf and pf are 
finite in the mixed phase - at low as well as high den- 
sities (or equivalcntly, bare chemical potentials). Hence 
the mixed phase in the weak coupling limit is described 
by a two-component Luttinger liquid at all densities. 

In the strong coupling limit, the behavior is markedly 
different at low and high densities. At low densities 
(jif + Tig < 1) the stiffnesses are finite in the mixed 
phase, impying a Luttinger liquid character of the ground 
state. At large pb/Fi the system is always at the com- 
mensurate filling, riF + ns = 1. In the pure phases, the 
stifnesses pf and pf vanish with the corresponding den- 
sities, ns and uf- Also, pf vanishes in the B-B phase in 
the thermodynamic limit since the bosons arc in the Mott 
insulating phase at this C/fl6|]. In the mixed phase, the 
individual stiffnesses have small, but finite, values . 
This suggests the existence of a gapless mode, in contrast 
to previous studies jL4|. The stiffness data and density 
correlations (shown below) support the mapping to an 
equivalent XXZ chain, which at U = V further reduces 
to the gapless XY chain. The gapless mode consists of 
the k = spin- wave created by the superoperator Q and 
corresponds to the "pairing superfluidity" in the BFHM 
explored in detail in Ref. [n]. 

To further characterize the nature of the mixed phase, 
we have studied the static charge susceptibility [Fig. |J] 
for the two components. Once again we find the results 
to be markedly different in the strong and weak cou- 
pling regimes. Fig. |](a) shows the results of increasing 
inter-species interaction on the static charge susceptibil- 
ity in the weak coupling regime. In the absence of boson- 
fermion interaction, xf(q) shows the familiar peak at 
twice the Fermi momentum {2k f)- As the inter-specics 
repulsion is increased, the peak at 2kp decreases and 
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FIG. 4: (a) The static charge susceptibility for the fermions, 
Xf(?), and bosons, Xb(q), with varying inter-species interac- 
tion strength, V. The peak at the Fermi wave vector is almost 
completely destroyed for V = 4. The densities of the two 
species are tib ~ 0.17 and uf ~ 0.24. The open (filled) sym- 
bols and dashed (solid) lines represent the bosonic (fermionic) 
data, (b) The suseptibilities in the strong coupling limit at 
commensurate filling, ns ~ 0.67 and rip ~ 0.33. 

more weight is transferee! to lower momenta. Finally at 
large values of V, the peak at the Fermi momentum is al- 
most completely washed out implying that in this regime 
the fermions no longer retain an independent character. 
On the other hand, in the strong coupling limit, at com- 
mensurate filling (tib + np = 1), the results are dia- 
metrically opposite. The boson density correlation de- 
velops a peak at twice the bosonic "fermi momentum" , 
2ks = 2tttib- In this limit of large U, the bosons behave 
effectively as hard-core bosons whose diagonal correla- 
tions are identical to those of spinless fermions. 

To conclude, we have studied the ground state phases 
of the BFHM on a ID optical lattice. The analytical 
expressions (J2J) , ^) for the boson-demixing boundary, 
constructed assuming a continuous transition, arc in an 
excellent agreement with the numerical results. For the 
symmetrical interactions, ts = ip, Ub = Vf, the F-F 
phase becomes locally unstable at [Af = fJ-B + h c (n, V) > 
Ub [Eq. (||)], any phase along the diagonal ub = Hf is 
gapless to boson-fermion conversion, while the fermion 
demixing transition happens strictly for (is > fJ-F- This 
establishes the stability of the mixed phase at all densi- 
ties and couplings. For strongly interacting bosons and 
fermions, the extent of the mixed phase gets progressively 
narrower at low densities. On the other hand, for weakly 
interacting particles, the homogeneous phase is stabilized 



over a wide range of chemical potentials. Investigation of 
the properties of the B-F phase at U = V reveals that at 
weak interactions, the mixed phase is a two-component 
Luttinger liquid at all densities. For strong interactions, 
this statement is valid at low densities, but breaks down 
at the commensurate total density, rip + ns = 1. Nu- 
merical data suggest a gapless mode for the ground state 
and strong coupling mapping to the XXZ model predict 
that at U — V this gapless mode consists of a k = spin- 
wave. In the BFHM, this corresponds to fluctuations of 
the mixture composition plj. We have also explored the 
effects of boson-fermion interaction on the density corre- 
lations of the two components and find that the results 
are dramatically different in the weak and strong cou- 
pling limits. These signatures should be readily observ- 
able in experiments, and, if confirmed, would indicate the 
emergence of novel many-body phases in these systems. 
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